function bits=golomb_data(x,th)
T=2^th;
bits = 0;
 for t = 1:length(x(:))
    bits = bits + floor(x(t)/T)+1+th;
 end